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Abstract 



In this paper we present the 1-loop perturbative computation of the renormaliza- 
tion constants and mixing coefficients of the lattice quark operators of rank three 
whose hadronic elements enter in the determination of the second moment of Deep 
^ ' Inelastic Scattering (DIS) structure functions. 

We have employed in our calculations the nearest-neighbor improved "clover-leaf" 
lattice QCD action. The interest of using this action in Monte Carlo simulations 
lies in the fact that all terms which in the continuum limit are effectively of order 
a (a being the lattice spacing) have been demonstrated to be absent from on-shell 
hadronic lattice matrix elements. We have limited our computations to the quenched 
case, in which quark operators do not mix with gluon operators. 

We have studied the transformation properties under the hypercubic group of the 
operators up to the rank five (which are related to moments up to the fourth of DIS 
structure functions), and we discuss the choice of the operators considered in this 
paper together with the feasibility of lattice computations for operators of higher 
ranks. 

To perform the huge amount of calculations required for the evaluation of all the 
relevant Feynman diagrams, we have extensively used the symbolic manipulation 
languages Schoonschip and Form. 
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1 Introduction 



This is the second of two papers addressed to the problem of computing in 
lattice QCD the 1-loop renormalization constants and mixing coefficients of 
the operators of rank two and three whose hadronic matrix elements respec- 
tively determine the first and second moment of Deep Inelastic Scattering 
(DIS) structure functions. In the first paper (hereafter referred to as [1]) we 
have reported the results for the case of the rank two operators. Here we will 
extend these results to the quark operators of rank three. 

Let us start with a brief introduction to light-cone physics and a summary of 
the improvement program in lattice gauge theories [2-5]. We refer the reader 
to [1] or to the review paper of Ref. [6] for more details. 



1.1 DIS and moments of the structure functions 



The quark operators whose matrix elements are related to the moments of the 
DIS structure functions can be written in the form [7-9] 



1 _ <-> 

)qS 



o 



= ^ 7{w Dm"- D m} (1 ± 75)^ 
1 — <-> , s \f 



tf-m} = ^ ^ D^-- - D„ n} (1 ± 7s)y i>, (1) 



where the A-^'s are flavor matrices. They are gauge invariant, have twist two 
and are traceless and symmetrized with respect to all Lorentz indices. S and 
NS superscripts refer to Singlet and Non Singlet flavor structures. 

In the unpolarized cross section the 75 contributions present in Eqs. (1) average 
to zero. The other contributions have hadronic matrix elements of the form 

{p\0^}. m \p) = A N (ii)p^ ■ ■ -p m + trace terms, (2) 
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where p is the subtraction point. They contain long distance contributions 
(non-perturbative physics) and are related to the iV-th moment of the DIS 
structure functions by the equations [10] 



J dx B XB^ r k (q 2 ,x B ) = C N +i(q 2 //x 2 )A N+ i(/x). (3) 
1.2 Improved lattice QCD 



The nearest-neighbor improved action (SW "clover-leaf" action) employed in 
these calculations is obtained by adding to the standard Wilson action [11] 
the Sheikholeslami-Wohlert [4] nearest-neighbor interaction termQ 

AS{ = -ig a 4 ~r ^n^VFn^n- (6) 



This term modifies the standard order g quark-gluon Wilson vertex 
(V) b p c (k,k') = -g (t A ) bc 



. a(k + k') p . a(k + k') P 
r sin + vy p cos 



(7) 



by adding to it the extra piece ( "improved vertex" ) 

(VX c (k, k') = -gJ-{t A ) bc cos a{k ~ k ' )p ]T a pX sin a(k - k') x , (8) 



where k and k! are the momenta of the incoming and the outgoing quark 
respectively and p is the Lorentz index carried by the gluon. The quark and 
gluon propagators, as well as the interaction between the quark current and an 
even number of gluons, turn out to be unmodified. To the order to which we 
will perform our calculations we will not need the expression of higher order 
quark-gluon vertices. 



5 Here F n>fiv is not the naive lattice "plaquette" 
1 



2igoa z 



{U n ,fiv U^ pu ), U n ^ v — Un^pUn+^^Un^^U^v, (4) 



but rather the average of the four plaquettes lying in the plane fiv, stemming from 
the point n: 

F n ,fj,u = ^ X/ Pri,iiu = g. a2 X] (U n ,fiu — U^ tflu ). (5) 

fiv=± °0 fiv=± 
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Besides adding the term (6) to the action, for consistency [5] one also has to 
rotate all the quark fields appearing in Green functions according to the rule 
(see Appendix A for the definition of the lattice covariant derivatives) 

^— ^{i+^py (9) 

As it has already been noticed in [1], the transformations (9) are the sources 
of very many algebraic complications in perturbative calculations. 



2 Structure of the operators 



In this paper we want to study the lattice renormalization of the rank three 
quark operator 

OU = ^7{^,S T} ^. (10) 



Its hadronic matrix elements are related to the second moment of the re- 
distribution of quarks inside the hadron. 0^ VT ^ is the (parity conserving piece 
of the) second operator in the list given in Eqs. (1). 

We will restrict ourselves to the quenched approximation, and consequently 
we will not have to consider the mixing between the quark operator (10) and 



the gluon operator J2 P Tr 



. Actually both the experimental and 



the numerical determinations of the x^-distribution of gluons are much harder 
to obtain than the corresponding quantities for quarks. 



2.1 General considerations 

The explicit expression of the "improved" quark operator, obtained by taking 
into account the rotations (9) on the fermion fields, is given by 

( U}) ! = 

^7{ M DuD t } i> - ?hv DuD T } $ - i>l{» DuD T } i> + ^7{p DuD T } $ 
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ar 
16 



DuD T }P DuD T }P l/> 



+^ p 7{/J DuD T} 1p + i) P 7{m 3„Sr} i)~i)p 1{» DuDr} j> 



a 2 r 2 



i> p 7{ M DuD T }p i> - i> P 7{^ DuD T }P $ 



-i) p 7{M D v D T }p i> + i> P 7{M DuD T }p 4> 



As can be seen, the three-index case introduces a new feature in the expression 
of improved operators, namely the appearance of the cross- derivative term, 

ip D fl D u if), which has to be defined according to the formula 



■ip n D^D U Vv, 



1 



4a 2 



;i2) 



The more usual terms ip D fl D u ifj, if) D^D U ip and ip D^D V ip have respectively 
the expressions 



(13) 



Aa 2 



(14) 



l\) n D^D U %p n 



Aa 2 



4>n. (15) 



The requirement that fixes the exact form of the cross-derivative term is that 

the right-derivative, D u , and the corresponding left-derivative, D u , should be 
obtained one from the other by an integration by parts, that amounts on the 
lattice to a shift of one site in the v (in this example) direction. 

Using the nearest-neighbor improved action, we want to compute to 1-loop in 
the chiral limit the matrix elements 



where \p) is a one-quark state of momentum p and vanishing (renormalized) 
mass. The renormalization constants and mixing coefficients of O q can be 
extracted from the knowledge of the amplitudes M M ^ T (p). 

2.2 Operator mixing 

Lorentz indices in Eq. (16) have to be chosen carefully by looking at the trans- 
formation properties of the resulting operators with respect to the hypercubic 
group, H(4) [12,13]. This group, which consists of the discrete proper rota- 
tions (no parity reflections^) in four dimensions, is the remnant of the proper 
euclidean Lorentz group on a discretized space-time. 

In moving from a continuum four-dimensional relativistic theory to its lattice 
version, the proper (euclidean) Lorentz invariance 0(4) is broken to the invari- 
ance under the proper hypercubic group H(4). For this reason representations 
that under the 0(4) group are irreducible become in general reducible under 
H(4) when the corresponding operators are written on the lattice. Therefore, 
some special care must be exerted to avoid unwanted mixings. 

The most dangerous situation occurs if one considers the operator O q with 
\i — v — t, because it can mix (and indeed it does) with the lower- dimensional 
operator 



with a power divergent (~ l/o 2 ) mixing coefficient. Any simulation carried 
out with an operator with three equal indices will require a delicate non- 
perturbative subtraction of the lower-dimensional operator. This ends up in 
very large statistical errors in the results of these simulations. 

6 In the following we will restrict ourselves to the transformation properties of the 
operators under the proper hypercubic group, as we will only consider unpolarized 
structure functions, in which case the 75 contributions drop out from Eqs. (1). 



M^ip) = {p\Ol T \p), 



(16) 



(17) 
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Another possibility is to take |i / 1/ / /i). In this way one selects an 
irreducible representation of H(4), and no mixing is anymore allowed. We 
will take as a typical component of the multiplet of operators transforming 
according to this representation the operator 0{ 123 }- The drawback of choosing 
H 7^ v 7^ r(^ /i) in a Monte Carlo simulation is that two components of 
the momentum of the hadron are bound to be different from zero and from 
each other. This can lead to significant systematical errors coming from the 
increased sensitivity to the granularity of the lattice. 

A third possible choice (and the best one) is fi = v ^ r. In this case only 
one component of the hadron momentum has to be taken different from zero. 
Martinelli and Sachrajda have shown that power divergent subtractions can 
be avoided by using the particular combination [14,15] 

°DIS = °{411} ~ -j^°\^2\ + °{433})> ( 18 ) 

which turns out to belong to an irreducible representation of H(4). Using the 
standard Wilson action, they estimated the renormalization constant of this 
operator in the hypothesis of tadpole dominance of lattice perturbation theory. 

Unfortunately, the fact that the operator (18) belongs to an irreducible rep- 
resentation of H(4) is not enough for it to be multiplicatively renormalizable, 
because in the decomposition of a rank three tensor into irreducible represen- 
tations, the representation to which the operator (18) belongs appears more 
than once. As a consequence, mixing can occur among operators belonging to 
these equivalent representations. 

In the next subsection we will present the complete solution to the problem of 
decomposing in irreducible H(4) representations the particular tensor products 
of representations that arise when one considers the (parity conserving pieces 
of the) DIS quark operators of rank 3, 4 and 5 appearing in Eqs. (1). 

2. 3 Lattice transformation properties of DIS quark operators 

The proper hypercubic group H(4) has 192 elements. Its 13 irreducible rep- 
resentations are listed in Table 1, where an extra subscript is used to label 
different representations with the same dimension. For example, the represen- 
tation labelled as li is the identity, while the representation labelled as 1 2 is 
the totally antisymmetric one. In Table 1 we also give the corresponding 0(4) 
representation in SU(2) <g> SU(2) notation, when it exists. We recall that the 
representations 3s, 36 and 4 2 can be obtained as tensor product of 1 2 with 33, 
3 4 and 4i respectively. With an abuse of nota tion, in the literature they are 
often denoted respectively as (1,0), (0, 1) and (|, \). 
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Table 1 

The irreducible representations of the hypercubic group (the corresponding 0(4) 
representations are also indicated) 

li 1 2 2 3i 3 2 3 3 3 4 3 5 3 6 4i 4 2 6 8 

(Q,o) - - - - (i,Q) (Q,i) - - (ih) - - - 

Table 2 

Multiplet structure for 4± (g> 4i <g) 4i 



typical index pattern 


multiplets 


123 


4i + 4 2 + 8 + 8 


112 


3 ■ (4i + 8) 


111 


4i 



We are interested in the decomposition of tensor products such as 4i<E>4i<8> - • - ® 
4i to which all the operators listed in Eqs. (1) belong. According to Table 1, 
the irreducible 4i representation of H(4) corresponds to the irreducible vector 
representation (\,\) of 0(4). 

The case of the rank three operator is easily worked out remembering the 
corresponding 0(4) decomposition formula 

A In A lx A In , A In n A 3, n ,3 L .3 3, 
<2' 2 } ( 2' 2 } ® ( 2' 2 } = 4 • ( 2' 2 } + 2 ' ( 2' 2 } + 2 ' ( 2' 2 } + ( 2' 2 } - (19) 

One finds for H(4) 

4i ® 4i ® 4i = 5 • 4i + 4 2 + 5 • 8. (20) 

These representations can be organized according to the patterns and number 
of indices having the same value, as shown in Table 2. In the first column we 
give a typical representative for each type of index pattern. 

It is straightforward although tedious to extend the above results to the more 
complicated cases of operators of rank four and five. The H(4) decomposition 
of the 256 components of the lattice rank four operator is given by 

4i <g) 4i ® 4i ® 4i = 5 • li + 1 2 + 5 • 2 + 10 • 3 X + 6 • 3 2 

+ 10 • 3 3 + 10 • 3 4 + 6 • 3 5 + 6 • 3 6 + 16 • 6. (21) 

The decomposition results in 75 multiplets belonging to ten different types of 
representations, which again can be organized according to the patterns and 
number of indices having the same value, as shown in Table 3. 
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Table 3 

Multiplet structure for 4i <g> 4i ® 4i <8> 4i 



typical index pattern 


multiplets 


1234 


li + 1 2 + 2 + 2 + 3i + 3i + 3i + 3 2 + 3 2 + 3 2 


1123 


6 • (3 3 + 3 4 + 3 5 + 3 6 + 6 + 6) 


1122 


3-(li + 2 + 3! + 3i + 3 2 ) 


1112 


4 • (3 3 + 3 4 + 6) 


1111 


li + 3i 


Table 4 

Multiplet structure for 4i <g> 4i ® 4i <8> 4i ® 4i 


typical index pattern 


multiplets 


11234 


10 • (4i + 4 2 + 8 + 8) 


11223 


15 • (4i + 4 2 + 8 + 8) 


11222 


10 • (4i + 8) 


11123 


10 • (4i + 4 2 + 8 + 8) 


11112 


5 ■ (4i + 8) 


11111 


4i 



Finally, the 1024 components of the lattice rank five operator decompose in 
171 multiplets of only three different kinds, 

4i <g> 4i ® 4i ® 4! <g> 4i = 51 • 4i + 35 • 4 2 + 85 • 8, (22) 
and their classification is reported in Table 4. 

2.4 Choice of the rank three operators 

It is easily seen that the operator (18) belongs to the eight-dimensional repre- 
sentation of H(4). Furthermore, we see from Eq. (20) that when one decom- 
poses O q in H(4) multiplets there arise five eight-dimensional representa- 
tions. The crucial point here is that the hypercubic group H(4) possesses only 
one single eight-dimensional irreducible representation (see Table 1). There- 
fore all the five 8's that come out in the decomposition of the rank three 
operators into H(4) multiplets can in principle mix among themselves. This 
is at variance with what happens in the case of the rank two operator, where 
the operator 0? 4 n is multiplicatively renormalizable because it belongs to the 
irreducible six-dimensional representation, 6, which appears only once in the 
decomposition 4i ® 4i = li + 3i + 3 3 + 3 4 + 6. 
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A natural question to ask is what are the equivalent representations that can 
really mix among themselves in the decompositions (20), (21) and (22). For a 
general answer see Ref. [16]. For our purposes it is enough to note here that 
radiative corrections can not only induce any permutation of the values of the 
indices, but can also make pairs of equal indices flip their common value. As 
an example, from On 2 radiative corrections may generate O211, O442, O323 and 
the like, but not O123 nor O114. Similarly Om can mix with O133, O414 and so 
on. But from O123 only O213 and permutations can arise, and not operators 
like O122 or 44 i. 

For the case of the operator (18) we thus expect mixing only among the three 
eight-dimensional representations that can be constructed by suitably arrang- 
ing the operator components that have two and only two equal indicesQ. To 
deal with the mixing of these three copies of the eight-dimensional represen- 
tation it is more convenient to consider the following three eight-dimensional 
representations^: 

Oln ~ 2^422 + ^433)5 ^422 — 2^411 + ^433)5 

^311 — 2^322 + ^344)5 ^322 — T^Pzil + ^344)) 



»') 01,-^(0^ + 0^), 



in) Of u -±(Ol 2i + Ol u ), ... (23) 



where neither symmetrization nor antisymmetrization of the indices is to be 
understood. The point here is that the operator (18), which is the symmetric 
combination of the three kinds of operators in Eqs. (23), actually mixes with 
the two other independent linear combinations of them that do not possess any 
special symmetry. In other words, the interaction may change the symmetry 
properties of these operators. 

In the continuum the transformation properties for the rank three operators 
are much simpler. In fact, the symmetric combination (18) belongs to the 

7 No mixing can occur with the operator 0{ 41 p because as we said it belongs to 

the representation 6, which does not appear in the decomposition (20). 

8 Apparently the number of components of each multiplet is 12 instead of 8. But 
in fact there are 4 constraints, because the sum of the four components of each line 
in (23) is zero. 
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sixteen-dimensional (§, §) irreducible representation of 0(4), which only ap- 
pears once in the decomposition (19), and consequently the operator (18) 
is multiplicatively renormalizable. It follows that in the continuum, unlike 
what happens on the lattice, the two four- component multiplets of operators 
(^{m}' ^{222}' • • •) an d (^{123}' ^{234}' • • •)> which are part of this representa- 
tion, have the same renormalization constants as the operator (18). We report 
in Table 5 the value of this common renormalization constant in the MS 
scheme. Vice versa no mixing occurs between the operator (18) and the two 
other multiplets of operators with which it mixes on the lattice, because they 
can be arranged in the two irreducible 0(4) representations (|, |) and (|, |). 

As for the operator 0{i23}, on the lattice it belongs to the 42 representation. 
This representation appears only once in the product 4i <S> 4i <S> 4i, implying 
that the operator 0{i23} is multiplicatively renormalizable. 

In this work we have computed the renormalization constant of O/123} and 
the full mixing matrix of the operators (23). We have also verified that, as 
expected, on the lattice the traceless operator 0{m} is not multiplicatively 
renormalizable, and that it indeed mixes with the three operators belonging 
to the 4i representations with two equal indices, as indicated in the second 
row of Table 2. 

2. 5 On the feasibility of higher moments computations 

We conclude this section with a few observations on the operators of rank four 
and five. 

In the case of the rank four operators, we see from Table 3 that there are 
two operators that are multiplicatively renormalizable. As in the rank three 
case, they have all indices different. One of these representations is the totally 
symmetric li representation, and the other is the completely antisymmetric 
I2 representation. The I2 representation is unique in the decomposition (21), 
while the li appears five times, grouped in three different patterns. But even 
in this case no mixing can occur, because if one starts with the pattern of 
indices 1234, radiative corrections can never lead to the other two kinds of 
patterns which the other four li representations possess, namely 1122 and 
1111. For any other choice of indices mixing will be unavoidable, because all 
the other representations are present more than once in the decomposition 
(21), and at the same time they have patterns of indices that can be turned 
one into the other by the interaction. 

The use of the rank four operator with all indices different in Monte Carlo 
simulations would require three components of the momentum of the hadron 
to be different from zero, and different among themselves. 
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For the rank five operators, the situation is of course even worse. In fact there 
are no representations that are present only once in the decomposition (22). 
Furthermore, we see from Table 4 that within a single pattern of indices (i.e. 
within a single row) a representation is never present only once, except for 
11111, which is a 4i. In this case, however, radiative corrections can induce it 
to mix with the ten 4i representations with index pattern 11222 and with the 
fifteen 4i representations with index pattern 11223. Therefore for rank five 
operators mixing is unavoidable. 



3 Renormalization conditions 



In the tree approximation the amputated non-vanishing matrix elements of 
the operators 0^ UT that we consider in this work are given by 



{p\Ol m (a)\p) 



tree 
amp 



1 



;7{iP2P3} 



(24) 



1 , „„ „„ s , s , , tree 1 / 2 1 / 2 2\ 

^ = -2^-2(74^+74%) 



(p\Ol 11 --(Ol 22 + Ol 33 )(a)\p) 



(p\Ol n -\{Ol i2 + Ot^{a)\p) ^ 



(P\0 q iu - 2^224 + O q 33A )(a)\p) 



amp 

tree 
amp 



~2 (71P1P4 - ^(72^2^4 + IsPsPi) 

~2 (7lPlP4 - ^{l2P2P4 + 12,PZP A) 



The renormalization conditions that connect the bare operators on the lattice 
to the finite operators renormalized at a scale \x must be the same as the 
ones imposed in the continuum. We define finite renormalized operators by 
requiring their two quark amputated matrix elements computed at p 2 = fi 2 
to be equal to the corresponding tree level amplitudes. The renormalization 
constant ^{123} for the multiplicatively renormalizable operator 0| 123 ^ is thus 
fixed in perturbation theory by simply imposing the renormalization condition: 



l°{123}(^)b) 



(p\°{123}( a )\P) 



tree 



(25) 



where ^{123} is given by 

<5{i23}(^) = Z{i23}MOf 123} (a). (26) 

The anomalous dimension of the multiplicatively renormalizable operator 0| 123 
is 25/3. For /ia <S 1, the renormalization constant of Or 123 | can be written in 
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the form 



'{123} M = 1 - log //a + 5 {123} 



(27) 



where a$ = 5 , o/^ 7r anc l Cp = (N% — 1)/2N C is the quadratic Casimir of the 
fundamental representation of SU(iV c ); in QCD N c = 3 and CV = 4/3. The 
constant -8(123} is all we have to know to fix Zs^ii,} and is given in Sect. 5 for 
different values of the Wilson parameter r. 

For the case in which two of the indices are equal, we see that two of the three 
possible combinations give the same tree level result (see Eqs. (24)). Therefore 
we can reduce the whole problem of mixing to the mixing of the two operators 

Oa = Ol u -^(Ol 22 + Ol 33 ) (28) 

and 

B = OLi + 0\ u - \ (Of 42 + Of 24 + O343 + O q 3U ). (29) 

If we set 

Oa = Z aa Oa + Z AB B 
B = Z BA A + Z BB B , (30) 



the mixing coefficients will be extracted from the 1-loop calculations according 
to the formulae 



(p\o A \ P ) 



1— loop 



1 



1 



1— loop 



+ 2Z AB (7lPlP4 - ^(72^2^4 + 13P3P4 

(P\0 B \P) 

+ 2Z BB (llPm - ^(l2P2P4 + 73P3P4))) • 



2 {^b\ [idv - \{liPl + 74P3) 



(31) 



Taking into account the form of the matrix of the anomalous dimensions of 
A and Ob (see Table 11), we can write 



(13 \ 
— log fia + B AA J 
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Z AB (na) = —t-C f (2 log im + B AB ) 

47T 

Z BA (fia) = ~-—C F (4 log \ia + B BA ) 

a$ ( 19 \ 
Z BB (/j,a) = 1 - ~^C F I— log/xa + B BB \ , 



(32) 



where the values of the constants B can be found in Sect. 5. 

For DIS the physically relevant operator is the symmetric combination O^is 
(Eq. (18)), which is simply given by 

O q ms = \(Oa + O b ) . (33) 



For it we have 



- {z AA + Z BA ) A + - (Z AB + Z BB ) B . 



(34) 



4 Some aspects of the computation 



The computation of the renormalization constant of the operator 0{i23} has 
been crosschecked by comparing the results obtained from two completely 
independent codes, written to automatically carry out all the extremely com- 
plicated algebraic manipulations necessary in the calculation [17]. In fact, 
contrary to the case of the rank two operator, no complete hand checks are 
possible hereQ: the operator of rank three contains one more covariant deriva- 
tive, and all propagators, vertices and operators have to be expanded to second 
order in a. This results in impossibly long algebraic expressions. In the case of 
the operator with two equal indices we have checked by hand the results only 
for the simplest Feynman diagrams. 

We started our computations by employing a Schoonschip code [18] which was 
obtained as an evolution of the one we used in Ref. [1]. At a later stage we 
decided to use the more flexible Form language, and we developed a completely 
new set of routines to deal with Dirac 7-matrices. 

9 We have completed the calculations by hand only for the non-improved diagrams; 
the calculations of the improved ones is hard drudgery. 
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The most innovative feature of the new Form program is a routine that reduces 
the products of 7-matrices in a much more efficient way than before. As a 
byproduct, the treatment of even products of sinesp ] can be postponed, and 
performed after the completion of the Dirac algebra; in this way, we happily 
no longer need to evaluate the daring expressions involving six and eight sines 
that gave rise to a plethora of terms (see Tables C.2 and C.3 in Ref. [1]) in 
the old Schoonschip programs. One discovers that what is needed in this new 
way of proceeding is at most the peculiar combination of six sines of the kind 

l(fii,fi 2 , fi 3 , ^4) = / d 4 k /(cos k, ^2 sin2 k\) sin fc w sin k^ sin k^ sin 3 k^ = 
J x 

/ d*k / (cos k, sin2 &a) sin 2 k^ sin 4 k^ 6 MljU2 5 
J x 

+ / d 4 k / (cos k, ^2 sin 2 k\) sin 2 sin 4 k^ 2 5 
J x 

+ d 4 k f (cos k,Y, 

sin 2 k\) sin sin 2 k^ 2 5 

J x 
-2- Jd A k / (cos k, sin 2 kx) sin 6 fc M1 5 M1(U2M3M4 , (35) 



where Sn lfl2ll:ifl4 is non-zero only if all the indices are equal. The structure and 
the coefficients of Eq. (35) are the same as those that appear in the similar 
formula for the case of four sines (see Table C.I in Ref. [1]), with a sine raised 
to the third power here playing the role of a plain sine there. 

The Form code described above is the one that has been carefully optimized 
and used more extensively in our computations. Meanwhile, another Form 
code, completely independent from it, has been developed by one of us [19]. 
Summarizing, three sets of routines have grown up from our efforts, and we 
have checked for each code separately that the rather complicated Dirac alge- 
bra (we have products of up to seven gamma matrices) is carried out correctly. 
The analytical results obtained with the modified old Schoonschip program 
and with the two newly developed Form codes have been compared, and they 
turn out to be in perfect agreement. This we regard as a significative check of 
our calculations. 

The numerical integration has also been performed using various different 
Fortran programs, that all gave the same results within a 1% relative error. 
For the numerical integration, and for the treatment of the logarithmic diver- 
gences, we have used the same strategy as in Ref. [1] and once more we refer 
the reader to that paper for details. 

10 Odd number of sines integrates to zero. 
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A big problem we have encountered was the limitation on the working memory 
allowed by Schoonschip and Form on the different machines we have used. This 
is of primary relevance in the first stages of the computation, when all vertices, 
propagators etc. have to be expanded up to second order in the lattice spacing 
a. When dealing with the rank three case, one encounters several products of 
up to ten trigonometric functions, each one to be expanded up to order a 2 . 
Every product at first will thus give rise to something like 10 4 monomials. A 
large part of them does not contribute to the final expression, to the order 
in a we are interested in, and has to be killed at the earliest possible stage 
of the calculation. A careless way of programming may lead to an increase of 
CPU times of several orders of magnitude, or even to an abrupt stop of the 
elaboration for lack of memory. 

In going to the new Form programs, and thanks to the evolution of the various 
routines, the CPU time needed to perform the whole analytic calculation was 
drastically cut down. The CPU time required by the old Schoonschip code 
running on a Sun 3 workstation ranged from a few minutes up to more than 
four hours for the most complicated diagrams. The Form program was first 
running on VAX- VMS machines and then on a HP-UX 9000/735 machine. In 
its present form the computation of a generic diagram only takes from about 
ten seconds up to a few minutes. 

We conclude this section with a discussion on an important physical point 
that needs some clarification. This is the role of the terms a 2 p 2 that appear 
in what we may call the effective "Feynman propagator", P(ak,ap;q), of a 
diagram, i.e. the overall denominator that arises after the introduction of the 
Feynman parameters. Generally it will have the form 

P(ak, ap; q) = g(ak) + ap^h^ak) + a(l - a)a 2 p 2 + a 2 p ll p v f liu (ak) 9 ,(36) 



with q > 2, where the functions h^(ak) and f^ u (ak) have the following small 
ak behavior: 

h^(ak) ° (ak) 3 
Uv(ak) afc ~° (ak) 2 . (37) 



We see that the terms a 2 p 2 enter either multiplied by a function of ak vanishing 
as ak — > (e.g. sin 2 ak, 1 — cos ak etc.) or by a constant independent of a (i.e. 
a(l — a)a 2 p 2 ). They should be treated differently according to whether the 
diagram has already the required p^p v tensor structure in front or whether 
this structure will emerge after expanding the Feynman denominator itself. 
For dimensional reasons in the first case the diagram will be logarithmically 
divergent. Terms such as ap^h^(ak) and a 2 p^p u f fMU (ak) are to be set to zero, 
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Table 5 

Values of the constant B 3 3 in the MS regularization scheme 
[222} 

sails -62/9 
vertex 4/9 

self-energy -1 
total - 67/9 

while the term a(l — a)a 2 p 2 must be kept in the denominator and will act as 
a regulator for the lattice infrared singularity at a k = 0. In the limit ap <C 1 
it will give rise to the expected logarithmic behavior, log ap. In the second 
case the integral is finite because a function of ak, vanishing as (ak) 2 in the 
limit ak — > 0, will dimensionally replace in the numerator the two missing 
momenta present in the previous case. In this situation the contribution of the 
diagram to the sought tensor structure is obtained by expanding the Feynman 
propagator in powers of a 2 p ll p v j '^(ak). Each term in this expansion is finite, 
thanks to the smoothing out provided by the function /^(a/c), and one is 
allowed to set to zero the a(l — a)a 2 p 2 term. Since we are interested in the 
terms that have at most two powers of p in front, we must look in the expansion 
only to the terms that are linear in a 2 p^p u f^ u (ak). 

Similarly P(ak, ap; q) must be expanded in powers of ap^h^ak) up to second 
order. Linear and quadratic terms must be kept in the expansion to combine 
with the possible p-dependence present in the numerator and reconstruct the 
required p^Pv tensor structure. The coefficients of this expansion are finite 
loop momentum integrals in which all a 2 p 2 terms can be safely set to zero. 



5 Results 

The results of our calculations are summarized in the Tables 5-12 presented in 
this section. The contributions coming from the standard Wilson action (that 
is, the non-improved results), those coming from the terms naively of order 
a in the improvement procedure and those coming from the terms naively 
of order a 2 are shown separately. Contributions from different classes of dia- 
grams, according to the classification given in Appendix B, are also presented 
separately. The accuracy of the numbers presented here depends only on the 
accuracy of the numerical integration and it is better than 1%. 

In Table 5 we give the continuum MS renormalization constants of the oper- 
ators 0{i23}, Oa + Ob and 0{m}. In Table 6 the results for the lattice renor- 
malization constants of the multiplicatively renormalizable operator 0{i23} 
are reported. Our numbers for the Wilson case agree, within errors, with 



r>MS 

{ -2'2> 
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Table 6 

Values of -B{i23} 





r 


Wilson 


0(a) impr. 


0(a 2 ) impr. 


total 




0.2 


-4.172 


-1.611 


-0.022 


-5.805 




0.4 


-5.368 


-3.683 


-0.106 


-9.157 


SAILS 


0.6 


-6.309 


-5.033 


-0.271 


-11.612 




0.8 


-7.032 


-5.838 


-0.520 


-13.389 




1.0 


-7.607 


-6.313 


-0.854 


-14.773 




0.2 


0.387 


-0.060 


0.012 


0.339 




0.4 


0.574 


-0.215 


0.064 


0.423 


VERTEX 


0.6 


0.804 


-0.398 


0.162 


0.568 




0.8 


1.021 


-0.572 


0.307 


0.756 




1.0 


1.216 


-0.731 


0.500 


0.984 




0.2 


-6.102 


-0.611 





-6.713 




0.4 


-4.326 


-1.513 





-5.839 


\ SELF-ENERGY 


0.6 


-2.762 


-2.311 





-5.073 




0.8 


-1.465 


-3.013 





-4.479 




1.0 


-0.381 


-3.646 





-4.027 




0.2 


-12.233 





0.009 


-12.224 


operator TADPOLE 


0.4 


-12.233 





0.037 


-12.196 


+ 


0.6 


-12.233 





0.084 


-12.149 


\ leg TADPOLE 


0.8 


-12.233 





0.150 


-12.083 




1.0 


-12.233 





0.234 


-11.999 




0.2 


-22.121 


-2.282 


-0.000 


-24.403 




0.4 


-21.353 


-5.411 


-0.005 


-26.769 


TOTAL 


0.6 


-20.500 


-7.741 


-0.025 


-28.266 




0.8 


-19.709 


-9.424 


-0.063 


-29.195 




1.0 


-19.005 


-10.690 


-0.120 


-29.815 



those of Ref. [20]. In Tables 7-10 we give, for r = 1 only, the values of 
the four lattice B constants relevant for the mixing between the operators 

OA = OZu-KQJb + QSss) andQ B = 0?41+0? 14 -K 242 + 2 9 24 + Ol43 + 3 9 34)- 

In Table 11 for completeness we report the mixing matrix of the anomalous 
dimensions of these operators. In Table 12 a summary of the mixing matrices 
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Table 7 



Values of Baa for r = 1 





Wilson 


0(a) impr. 


0(a 2 ) impr. 


total 


SAILS 


-2.353 


-3.568 


-5.582 


-11.502 


VERTEX 


0.870 


-0.258 


0.384 


0.996 


\ SELF-ENERGY 


-0.381 


-3.646 





-4.027 


operator TADPOLE 










+ 


-16.960 





4.104 


-12.856 


\ leg TADPOLE 










TOTAL 


-18.824 


-7.472 


-1.093 


-27.389 


Table 8 

Values of Bab for r = 1 




Wilson 


0(a) impr. 


0(a 2 ) impr. 


total 


SAILS 


-0.953 


-1.183 


2.391 


0.255 


VERTEX 


0.029 


-0.245 


0.007 


-0.208 


\ SELF-ENERGY 














operator TADPOLE 










+ 








-3.649 


-3.649 


\ leg TADPOLE 










TOTAL 


-0.924 


-1.428 


-1.251 


-3.603 


Table 9 

Values of Bba for r = 1 




Wilson 


0(a) impr. 


0(a 2 ) impr. 


total 


SAILS 


-3.108 


-4.731 


5.465 


-2.373 


VERTEX 


0.153 


-0.556 


-0.044 


-0.448 


\ SELF-ENERGY 














operator TADPOLE 










+ 








-8.982 


-8.982 


\ leg TADPOLE 










TOTAL 


-2.955 


-5.287 


-3.562 


-11.803 
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Table 10 

Values of Bbb for r = 1 





Wilson 


U{a) impr. 


(J(a ) impr. 


total 


SAILS 


-6.421 


-3.899 


-4.313 


-14.633 


VERTEX 


1.494 


-0.462 


0.725 


1.758 


1 OT71T T71 T71 TVTTT1 T* /~1 ~V7~ 

2 SEL1 1 -ENERGY 


-0.381 


-3.646 





-4.027 


operator TADPOLE 










+ 


-12.233 





10.596 


-1.637 


\ leg TADPOLE 










TOTAL 


-17.540 


-8.006 


7.009 


-18.538 



Table 11 

Mixing matrix of the anomalous dimensions of the operators Oa and Ob 



ANOMALOUS DIMENSIONS 



SAILS 



10/3 5/3 
10/3 5 



VERTEX 



-1 1/3 
2/3 -2/3 



i SELF-ENERGY 



2 
2 



TOTAL 



13/3 2 
4 19/3 



for the B constants both in continuum MS and on the lattice (where we give 
separately the results for the Wilson case and for the fully improved case) is 
given. 

Two observations are in order here. First of all we notice that for the physically 
relevant symmetric operator Op IS = \{Oa + Ob) the sum of the elements 
in each column of the lattice matrices are not equal, signalling that mixing 
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Table 12 

Mixing matrix for the B constants (the lattice values are for r = 1) 



Continuum MS 



Lattice 
(Wilson) 



Lattice 
(full improved) 



SAILS 





31 


31 




9 


18 




31 


93 


9 


18 



-2.353 -0.953 
-3.108 -6.421 



^ -11.502 0.255 \ 



-2.373 -14.633 



/ 



VERTEX 


1* 

9 


I) 


^0.870 




2 




0.153 


SELF-ENERGY 




:) ( 


-0.381 













0.996 -0.208 
-0.448 1.758 

-4.027 
-4.027 



TADPOLES 



V 



-16.960 
-12.233 



V 



-12.856 -3.649 
-8.982 -1.637 



TOTAL 




-18.824 -0.924 
-2.955 -17.540 



^ -27.389 -3.603 \ 



-11.803 -18.538 



/ 



occurs. In the continuum instead, as expected from the discussion of Sect. 2, 
it can be easily checked from the numbers given in Tables 11 and 12 that 
the combination O a + Ob is in fact multiplicatively renormalizable. A second 
observation that emerges by looking at the results of Table 12 is that, at least 
in this case, one should talk of "tadpole + sails dominance" rather than simply 
of "tadpole dominance" in lattice perturbation theory. 

We may summarize the results obtained in this paper in the following formulae. 
1) Standard Wilson action (Tables 6 and 12): 



^{123} 



25 , 1 
1 - —77-: log jia - 773^ ^{123} 



1 - 



6tt 2 /? 
0.422 



log jia + 



0.963 



•^{123} 



^{123}' 



(38) 
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(l q 



1 

+ 3 



1 - 



25 



log \ia 



6tt 2 /3 & ^ 2^ ft 



~^b( Baa + Bba 



O a 



25 1 



1 0.422 , 
3 1 -^- l0 §^ + 



6ir 2 ft 
1.103 



P 



A + 



2n 2 ft 
1 



0/ 



0.422 , 0.935 
1 — log jia + 



ft 



ft 



Ob- (39) 



2) Nearest-neighbor improved action (Tables 6 and 12) 



(0{123}Y 



25 1 



(0{123}Y 



0.422, 1.510 
1 — log fia + 



P 



P 



(OI123}) 1 ' 



(40) 



(Obis) 1 = 1 



25 
6^p 



lo£ 



(OaY 



+- 



1 - 



25 

6^(3 



lo g^ a - ^p{ B AB + B BB 



(Ob) 1 



0.422, 1.985 
1 — log fia + 



+ 



P 

0.422 



■ log jia + 





1.122' 
~P~ 



(OaY 



(ObY- 



(41) 



Obviously physical quantities should at the end be independent from the cho- 
sen regularization procedure. In particular the dependence from the subtrac- 
tion point, /i, must disappear from physical hadronic matrix elements. In fact, 
consistently to each order in perturbation theory, the log /ia terms get can- 
celed in the product between the Wilson coefficients and the matrix elements 
of the renormalized operators that are eigenvectors of the anomalous dimen- 
sion matrix. The net result is that effectively the Wilson coefficients must be 
taken at a momentum scale a" 1 and the operators turn out to be effectively 
renormalized by "reduced" renormalization constants Z, obtained from the 
full expressions by dropping all log \xa terms. 

Numerically for N c = 3 at the typical values ft = 2 N c / Qq = 6 and r = 1 and in 
the quenched approximation, one gets for the unmixed quark operator 0| 123 ^ 
from Eqs. (38) and (40) 
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0| 123} = 1.160 0\ m Wilson case 
(^{m}) 1 = 1-252 (O q 123 yY Improved case , 

where by the superscript we mean that the log //a term has been dropped 
from the expression of the 2Ts, as explained before. 

Similarly in the case of the operator 0^ IS one has from Eqs. (39) and (41) 

0£is= \ [1.184 0a + 1.156 B ] Wilson case 

(d q Dls Y = | [l.331 (0 A y + 1.187 (Ob) 1 ] Improved case . 

In the existing simulations corrections due to renormalization and mixing ef- 
fects were either calculated by means of tadpole dominance [14,15] or for 
simplicity neglected [21]. We see that even this latter approximation is not 
quite adequate because the effective renormalization constants and mixing co- 
efficients are not small. It should also be noted that the renormalization con- 
stants do not seem to decrease in magnitude with the order N of the moment. 
This fact is already evident in the Wilson case, where the renormalization 
constants of the rank three operators are larger than those of the rank two 
quark operators (see Ref. [1]), and it is confirmed by the complete improved 
computation. 



A Perturbative expansion of the vertex operator (O 



q Y 



In this appendix we present the perturbative expansion of the vertex operator 
{O q ) l , given by Eq. (11). We have used the definitions 



D n Ipn = 7T - [U„,v1pn+» ~ Un-^n-t] (A.l) 



2a 
1 

2a 



and the conventions 



A n „ = J^ 4 j {q + q » /2)n AM (A.2) 
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where the integrals are performed over the first Brillouin zone. Throughout 
this appendix external and loop momenta are expressed in lattice units. 



Since the full Fourier transform of the operator (0^ T ) r is very complicated, 
we give here only the form it effectively takes when inserted in the diagrams 
of Fig. B.2. Calling p the external incoming and outgoing momentum, and k 
the fermion loop momentum, one finds, separating the various contributions 
according to their naive order in a: 



a) tree level 



iar r d 4 k — . sin k u sin k v sin k T , .,. , . „s 

+ t / os^Wnr— (A - 5) 



(2*) 



a 2 r 2 r d 4 k _ „ sin k\ sin k x > sin k v sin k T , , , . 

+ / TTVa^ fc 2- 7a7 m 7a' ^(fc). 

8 J (27r) 4 aa' ° a a a 



b) order g 



(O^ T ) I (n = 0)| 



go r d 4 k — 



90 



2 J (2tt) 



£; + p sin A; 
cos(— — )„ A„(p - fc) 



£: + p ship,, i 
+ cos( — — ) T A T (p-k) ip{k) 



iag r f d 4 k 



(27T) 



cos 



fc + p sin /c r sin fc> 



a a 



A v {p - k) 



k + p sink x smp u 

+ cos(— — ) T A T (p - k) 

1 a a 



k + p ship,, sin p T 

+ cos (^— )a A x (p - k) 

2 a a 



r ( k + p sin fc T sin p x , M 

+7a7 m cos( — — )„ A,(p-fc) 

L 2 a a 

A; + p, sin sin pa 



+ cos 



+ cos( 



2 a a 
k + p sin fc„ sin fc T 



A T (p - k) 



a a 



A x (p-k) 



(A.6) 
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Fig. B.l. The graph that symbolically represents the 1-loop vertex correction to the 
insertion of the (O^) 1 operator. The operator insertion is indicated by a dot. The 
wavy line is a gluon. 



a 2 g r 2 



r d^k — 

W ,(p, § 



,k + p sin p\r sin k T sin k\ 
cost - )v A v [p - k) 



a a 



a 



k + p smp v smpxsmk x . , M 

+ cos(— — ) T A T {p - k) 

1 a a a 

,k + p sin k u sin k T sin ky . . 

+ cos(— — ) A A x {3> - k) 

l a a a 



,k + p sinp„sinp T sinpA' „ 

+ cos(— — ) x > A y {p - k) 

l a a a 



The above formula is given in the kinematical configuration in which the 
incoming gluon momentum lands on the incoming quark leg. If the gluon is 
attached to the outgoing quark leg, one must exchange p and k. Notice that the 
first lines in Eqs. (A. 5) and (A. 6) correspond to the non-improved expression 
of the operator. 

We do not give here the expression of the O(g^) terms because, besides being 
extremely complicated, they are not actually necessary for our computation. In 
fact we need them only when either the gluon or the quark legs are contracted 
to make a tadpole loop, and in this situation the tadpole directly factorizes 
out. 

The case in which two or more indices are identical is easily derived from the 
expressions above. 
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(a) Vertex 




(e) Self-energy 




(g) Leg tadpole 

Fig. B.2. The different types of graphs 
the matrix element {p\(O q ) l \p) . 



(b) Sail 




(d) Operator tadpole 




(f) Self-energy 




(h) Leg tadpole 



contributing to the 1-loop approximation 
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B Diagrams 

In what follows we show the 1-loop diagrams that we have calculated for each 
one of the operators discussed in this work. They are accompanied by a label 
that allows an easy connection with the Tables presented in Sect. 5. 

Actually, each graph in this appendix corresponds to several diagrams in the 
perturbative expansion. For example, the vertex correction of Fig. B.l symbol- 
ically represents the sum of 12 different diagrams. They come out remembering 
that for the quark-gluon interaction one can either take the standard Wilson 
vertex of Eq. (7) or the improved vertex of Eq. (8), and for the operator in- 
sertion either its unrotated form or the 0(a) or else the 0(a 2 ) corrections (see 
Eqs. (11) and (A.5)). 

We report in Fig. B.2 the types of graphs that contribute to the 1-loop cal- 
culation of the matrix element (p\(O q ) l \p), with the understanding that the 
drawings we show are representative of the structure of full sets of improved 
diagrams. 



C a-integration 

To reduce the computing time in the numerical evaluation of the integrals 
that involve fermion propagators, we have chosen to perform analytically the 
integration over the Feynman parameter a, and to leave for the numerical in- 
tegration only a four-dimensional expression. To this end after the subtraction 
of the standard logarithmic divergence we need to compute integrals of the 
form 



i 



F nm (f(k),g(k)^) = J 



a 



,71 



da 



[f(k) + ag(k)} 



(C.l) 



o 



In our calculations we explicitly have 



A 




[$>in 2 k x + 4r 2 C£ 



A A 





(C.2) 



27 



The functions F nm satisfy the recurrence relations 



of 



d_ 

dg 



F 77,777, 777- ' F 



n+1 m+1 



(C.3) 
(C.4) 



which make simpler their computation and allow for a check of the formulae 
given below. 

The formulae needed in this work are: 



F 2 : 

F 12 = 
F22 ; 
F 32 : 
F42 : 

F 3 : 
Flv 

F23 : 

F 3 3 : 
F43 : 
F 5 3 : 
F 6 3 : 
F 4 : 

F M = 



/•(<?+/) 

1 

9 



H l + 7. 



9 + f 



1 

7 

1 



9 Z 



/ 2 



2 
"1 

3 

9 + V 



g-(g + f) 



_ Z + 3 £ - 4^ log f 1 + 7) + 2 A n 
5 # 2 g 3 \ f g 2 -{g + f) 



2/ 2 • (9 + f? 
1 

V-{g + f? 



g £ 
1 

1 



3^ + 2/ 
2-(<7 + /) 5 



H 1 + l 



3/ , J 2 



J 3 



g 



9g + 8f 



3-2^ 6 7- 10 7 lo H 1 + 7j +/ 'v-(. + /) 2 



1 / J 2 



g 



g 



10 7 + 1 iH 1 + T 



-f 



4 iifl + 10/ 

2^ 3 • (g + /) 5 



£ 2 + 3gf + 3/ 2 
3/ 3 • (g + f) 3 

g + 3f 
6/ 2 • (g + /) 3 



(C.5) 
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1 

= J_ [lj ( 1+ 9\ 25g 3 + 52g 2 f + A2gp + 12f 





■ j >~\ 
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